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Abstract.  The  merging  of  an  intense  localized  vortex  with  an  extended  vortex  is  investigated  by 
means  of  an  experimental  analysis  performed  on  a Malmberg-Penning  trap  with  photocathode, 
and  numerical  simulations  with  a two-dimensional  (2D)  particle  in  cell  (PIC)  code.  The  study  is 
restricted  to  highly  nonlinear  conditions,  where  the  perturbative  approach  does  not  hold.  A very 
good  agreement  between  experimental  results  and  simulations  is  obtained.  It  is  found  that  the 
localized  vortex  is  firstly  wrapped  around  by  the  extended  vortex,  then  moves  towards  the  center 
of  the  system,  eventually  approaching  an  almost  stationary  state  of  rotation  characterized  by  the 
formation  of  vorticity  holes.  During  the  whole  evolution  the  extended  vortex  gains  energy  from  the 
field  surrounding  the  vortices,  while  the  energy  of  the  localized  vortex  remains  constant. 


INTRODUCTION 

Under  suitable  conditions,  the  dynamics  perpendicular  to  the  imposed  magnetic  field 
of  a pure  electron  plasma  confined  in  a Malmberg-Penning  trap  is  equivalent  to  that  of 
a 2D  incompressible,  inviscid  fluid  [1],  where  the  vorticity  C,  is  proportional  to  the  2D 
electron  density  n.  This  fact  allows  to  investigate  on  such  a system  nonlinear  processes 
occurring  in  2D  inviscid  flows;  in  particular  vortex  merger  and  filamentation,  which  play 
a dominant  role  in  the  evolution  of  turbulence  [2]. 

We  consider  here  the  process  of  interaction  between  an  extended  vortex  and  one 
very  localized  vortex  (almost  pointlike)  under  conditions,  which  have  now  become 
experimentally  accessible,  where  the  existing  perturbative  analysis  [3]  does  not  apply. 
The  experimental  investigation  of  this  process  is  performed  on  a Malmberg-Penning 
trap  (see  Fig.  1)  with  photocathode  [4],  which  allows  to  produce  well-defined  and  well- 
localized  vortices  of  almost  uniform  vorticity.  The  interpretation  of  the  process  is  based 
on  the  Hamiltonian  character  of  the  motion,  and  numerical  simulations  with  the  2D  (in 
real  space)  PIC  code  XOOPIC  [5]. 

Numerical  investigations  of  the  interaction  of  two  isolated  vortices  of  uniform  vortic- 
ity, were  reported,  e.g.,  in  Refs.  [6,  7,  8],  In  Ref.  [6],  in  particular,  it  was  pointed  out  the 
possibility  of  a merger  by  “entrainment”  of  the  larger- vorticity  region  within  the  smaller- 
vorticity  region.  In  Ref.  [8],  all  possible  regimes  of  interaction  were  investigated,  but 
assuming  identical  vorticity  for  the  two  vortices.  These  works  were  all  based  on  contour 
dynamics  simulations.  In  Ref.  [7],  based  on  the  use  of  a low-order  Hamiltonian  moment 
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FIGURE  1.  The  Malmberg-Penning  Trap  with  the  photocathode  [4],  The  desired  initial  2D  electron 
distribution  is  created  by  projecting  the  appropriate  light  image  onto  a cesium  antimonide  photocathode 
and  grounding  the  left  cylinder;  electrons  are  emitted  only  where  there  is  light,  and  they  stream  along  the 
magnetic  field  lines  into  the  central  confinement  region,  preserving  their  distribution. 


model  and  a pseudospectral  algorithm,  it  was  pointed  out  how  the  interaction  between 
two  unequal  vortices,  of  different  size  and  vorticity,  could  cause  one  of  the  two  vortices 
to  dominate.  In  the  language  of  Ref.  [7],  the  case  studied  here  corresponds  to  the  smaller 
vortex  being  victor. 

The  dynamics  of  the  system  is  described  by  the  divergence-free  velocity  field 
u(jt,y,f)  = e,  x Vt|/(x,y,/),  where  the  evolution  of  the  stream  function  \]/  is  determined 
by  the  conservation  of  the  vorticity  £ h e:  • V x u, 

dCM  + M = o,  CD 

where  [a,b\  = da/dy  ■ db/dx  - da/dx  ■ db/dy  is  the  Poisson  bracket,  the  coordinates  y,x 
are  conjugate  variables  with  respect  to  the  Hamiltonian  \j/,  and 


V2t|/  = c (2) 

The  functions  \j/  and  £ are  related  to  the  electric  potential  <|>  and  electron  density  n by 
the  relations  <j>  = (B/c)\\i  and  n = (B/4nec)C„  respectively,  where  B is  the  magnetic 
field  strength  and  the  other  notations  are  standard.  The  evolution  of  the  system  is 
determined  by  appropriate  initial  and  boundary  conditions  for  xj/.  In  particular,  the  free- 
slip  boundary  condition  \\i  = 0 on  the  circular  wall  corresponds  to  the  condition  <|>  = 0 
for  the  electrostatic  potential. 

We  refer  to  two  initially  circular  vortices  with  uniform  vorticities  <<  £p,  and  ini- 
tial radii  ae»  ap  (the  indices  e,p  denote  the  extended  and  the  pointlike  vortex,  respec- 
tively). The  experimentally  investigated  range  of  parameters  in  general  corresponds  to 
conditions,  in  which  the  initial  field  \|/  exhibits  a separatrix,  with  its  X-point  located  out- 
side the  two  vortices.  Reference  parameters  can  be  taken  Tp/Te  ss  10“ 1 , with  the  initial 
distance  d between  the  centers  of  the  two  vortices  in  the  range  1 + rp/r,,  < d/ae  < 2, 
with  re  p = n(^)Pa2e  p the  circulation  of  the  two  vortices.  A discussion  of  the  characteris- 
tics of  the  merging  process  when  the  above  parameters  are  varied  is  reported  in  Ref.  [9]. 

In  the  cases  under  consideration  the  analysis  which  can  be  canned  on,  following  Refs. 
[3,  10,  1 1],  has  a very  narrow  domain  of  validity  since  the  shape  of  the  extended  vortex 
is  soon  highly  distorted  from  its  initial  condition,  the  merging  being  very  fast  and  lasting 
few  rotation  periods. 
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FIGURE  2.  Time  evolution  of  the  interaction  of  a point  vortex  with  an  extended  vortex.  First  row: 
density  (experimental  results).  Second  row:  density  (PIC  simulation).  Third  row:  contour  plot  of  the 
corotating  potential  rjf  (PIC  simulation):  the  thick  line  represents  the  separatrix  crossing  the  extended 
vortex.  The  lengths  are  normalized  over  Rw.  From  left  to  right,  the  data  refer  to  t = 0.0,  0.5,  1.0,  and  4.0, 
respectively,  with  the  time  measured  in  units  of  the  period  of  rotation  of  the  extended  vortex,  4n/t^,. 


Fig.  2 in  the  first  row  shows  the  experimentally  detected  images  (the  intensity  is 
proportional  to  the  axially  integrated  electron  density)  at  different  times  of  the  evolution: 
from  the  initial  condition  to  the  final  almost  stationary  state.  The  corresponding  plots 
obtained  by  PIC  simulations  are  shown  in  the  second  row  of  the  same  figure.  The  number 
N of  macroparticles  used  in  the  simulations  is  < 105  with  a 256  x 256  square  grid.  The 
parameters  of  the  simulation  are  the  same  as  in  the  experimental  situation.  The  stability 
of  the  results  with  respect  to  variations  of  N and  grid  size  has  been  verified. 

The  agreement  between  experimental  data  and  simulation  shows  that  the  behavior  of 
the  system  is  indeed  2D  Eulerian,  even  in  this  highly  nonlinear  regime;  this  confirms  the 
potentiality  of  Malmberg-Penning  traps  for  fluid  dynamics  experiments. 

To  describe  the  evolution  of  the  system,  we  focus  on  the  behavior  of  the  stream 
function  in  a reference  frame  which  rotates  with  the  angular  velocity,  to,  of  the  localized 
vortex,  vjf  = vg  — o^x2  +y2) / 2.  The  initial  condition  of  the  system  appears  far  from  the 
equilibrium,  which  is  characterized  by  a vorticity  field  depending  on  the  space  variables 
only  through  t(7.  In  the  case  under  consideration,  we  chose  a value  of  to  corresponding 
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to  the  angular  velocity  of  the  localized  vortex,  0)  ~ Te/{2nd2).  For  high  enough  values 
of  the  ratio  Vp/Te  (verifying  the  conditions  mentioned  above),  a separatrix  of  the  tjf 
field  exists,  with  an  hyperbolic  point  located  between  the  centers  of  the  two  vortices,  as 
shown  in  the  first  plot  in  the  third  row  of  Fig.  2.  Moreover,  such  separatrix  crosses  the 
extended  vortex  so  that  a finite  part  of  it  lies  outside  the  separatrix:  this  part  gives  rise  to 
the  wrapping  process. 

The  region  enclosed  by  the  separatrix  around  the  localized  vortex  can  be  crossed  by 
few  fluid  particles  only,  and  forms  a zero  vorticity  region  which  surrounds  the  localized 
vortex  during  the  whole  evolution  of  the  system.  This  has  been  remarked  also  in  the 
analysis  of  the  motion  of  a localized  vortex  in  a space  dependent  vorticity  background 
[10],  The  presence  of  the  mentioned  hyperbolic  point  determines  the  shape  of  the 
extended  vortex  during  the  initial  phase  of  wrapping,  as  shown  in  the  second  plot,  a 
key  role  being  due  to  the  different  behavior  of  fluid  particles  in  the  neighborhood  of  that 
point.  As  wrapping  develops,  the  related  shift  of  the  center  of  vorticity  of  the  extended 
vortex  feeds  the  process  and  leads  to  the  coalescence  of  the  hyperbolic  point  with  the 
elliptic  point  associated  to  the  extended  vortex.  Then,  \j7  exhibits  a single  maximum  in 
the  region  of  interest,  as  shown  in  the  last  two  plots. 

The  second  phase  of  the  motion  is  characterized  by  formation  and  stretching  of 
filaments,  and  effective  reconnection  of  the  vorticity  field  with  the  related  formation 
and  stirring  of  holes,  which  occurs  due  to  the  Kelvin-Helmholtz  instability  arising  in  the 
spiral  channel  of  zero  vorticity  inside  the  extended  (constant  vorticity)  vortex  [7).  During 
this  process  the  \pr  field  does  not  show  significant  large  scale  evolution.  The  constraint  of 
vorticity  conservation  does  not  allow  the  system  to  reach  a state  where  £ is  a function  of 
\j7  only. 

Finally,  we  consider  the  behavior  of  the  energy  flows  during  the  merger.  The  energy 
density  is  given  by  £(*,  v,f)  = j|V\i/|2.  This  is  initially  concentrated  near  the  boundaries 
of  the  two  vortices,  and  it  has  a local  minimum  in  the  region  around  the  above  mentioned 
X-point.  The  total  energy  of  the  system,  £,  is  a constant  of  motion  and  reads 

E=l-fs\V\V\2dS  = '£E‘>  (3) 

where  the  integration  is  over  the  whole  space  within  the  circular  boundary,  and  the 
energies  £,,  with  i = e,p,b,  have  been  introduced,  representing  the  energy  associated 
to  the  extended  voitex,  the  localized  vortex,  and  the  region  surrounding  the  two  vortices 
(background),  respectively.  The  behavior  of  the  energy  amount  in  the  two  vortices  and  in 
the  background  is  shown  in  Fig.  3,  where  the  time  evolution  of  the  £,•’ s is  shown.  During 
the  whole  process  the  energy  of  the  pointlike  vortex  approximately  does  not  change, 
while  the  energy  of  the  extended  vortex  increases  (almost  linearly),  at  the  expense  of  the 
energy  of  the  background. 
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FIGURE  3.  Time  evolution  of  the  total  fluid  energy  of  the  two  vortices  and  the  background  (PIC 
simulation).  The  data  are  the  same  as  in  Fig.  2. 
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